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Abstract 

We present our recent results on the fragmentation by varying the mass asymmetry of the reaction 
between 0.2 and 0.7 at an incident energy of 250 MeV/nucleon. For the present study, the total 
mass of the system is kept constant (Atot = 152) and mass asymmetry of the reaction is defined by 
the asymmetry parameter (n = \ (At — Ap)/(Ax + Ap) |). The measured distributions are shown 
as a function of the total charge of all projectile fragments, Zb oun d- We see an interesting outcome 
for rise and fall in the production of intermediate mass fragments (IMFs) for large asymmetric 
colliding nuclei. This trend, however, is completely missing for large asymmetric nuclei. Therefore, 
experiments are needed to verify this prediction. 
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I. INTRODUCTION 



Heavy-ion collisions have always played a fascinating role in exploring various aspects 
of nuclear dynamics such as fusion-fission, multifragmentation and particle production. 
Multifragmentation, that is the emission of several intermediate mass fragments IMF's from 
a hot compound nucleus, is a phenomenon observed in nuclear reactions over a wide incident 
energy range. There has been considerable progress during recent years in the experimental 
studies. Experimental evidence for the statistical property of nuclear fragmentation has 



been given 



-|4|]. These quantities 



U-|3| and various new quantities have been measured 
include the mean multiplicity of intermediate mass fragments (< Njmf >), the average 
charge of the largest fragment (Z max ), the sum of all charges with Z > 2 etc. The 
quantity which is intimately related to the multifragmentation process is the multiplicity of 
intermediate mass fragments. Correlation between mean multiplicity of IMF's, < Njmf >, 
and the mass of the fragmenting system, whose measure is so called bound charge Z bound is 



an important aspect of multifragmentation that has been studied thoroughly by a number 
Q, Q, S]- They, however, didnot take asyn 



of groups |2|, |4|, |5[. They, however, didnot take asymmetry of the system into account which 
is very important to study the isospin effects jf], Q] ■ The asymmetry of the reaction can be 
defined by the parameter r\ = \ (A T — A P )/(A T + A P ) |; where A T and A P are the masses 
of target and projectile. The rj = corresponds to the symmetric reactions, whereas, 
non-zero value of rj define different asymmetry of the reaction. It is worth mentioning that 
the reaction dynamics in a symmetric reaction (rj = 0) can be quite different compared 
to asymmetric reaction {rj ^ 0) [8j. This is due to the deposition of excitation energy 
in the form of compressional energy and thermal energy in symmetric and asymmetric 
reactions, respectively. The multifragmentation is studied many times in the literature 
6|, |7| . Unfortunately, very little study is available for the mass asymmetry of the reaction 
in terms of multifragmentation. 

In recent years, it has become possible to do exclusive measurements of multifragmentation 



2 



process. This has been done with streamer chamber detectors, electronic detectors, and 

Att detectors. ALADiN 2|, |4j, |5| group has reported that the mean multiplicity of IMF's 

< Nimf > w as found to be same for all targets ranging from Beryllium to Lead and for 

r~l 
i i 

E/A ranging from 400 to 1000 MeV/nucleon. De Souza et al, |9[ observed a linear increase 
in the multifragmentation of IMF's for central collisions with incident energies varying 
between 35 and 110 MeV/nucleon. In 2009, Tsang et al, 10] reported a rise and fall in 



the production of IMF's. The maximal value of the IMF's shifts from nearly central to 
peripheral collisions with the increase in the incident energy. 

Theoretically, multifragmentation can be studied by statistical fl a S well a S dynamical 
models [12j]. The universal property of multifragmentation has been quite satisfactorily 



111 ] . On the other hand, dynamical 



described by the statistical multifragmentation models 
models are very useful for studying the reaction from the initial state to the final state 
where matter is fragmented and cold. In this paper, we will adress the most interesting 
dependence of the multiplicity of intermediate mass fragments (IMF's). This multiplicity is 
estimated in terms of the "bound" charge value. We have used Isospin-dependent quantum 
molecular (IQMD) model to study the effect of asymmetry of colliding nuclei on the 
mult if r agment at ion . 



The isospin-dependent quantum molecular dynamics 
charge states of nucleons, deltas and pions explicitly 



(IQMD) 13] model treats different 



14j , as inherited from the Vlasov- 



Uehling-Uhlenbeck (VUU) model |15|. The details about the elastic and inelastic cross 



sections for proton-proton and neutron-neutron collisions can be found in Refs. |13 
In this model, baryons are represented by Gaussian-shaped density distributions 



16|. 



fi(r,P,t) 



1 



TT 2 h 2 



2L P fi2 



(1) 



Nucleons are initialized in a sphere with radius R = 1.12 A 1 ^ 3 fm, in accordance with the 
liquid drop model. Each nucleon occupies a volume of h 3 so that phase space is uniformly 



filled. The initial momenta are randomly chosen between and Fermi momentum pp. 
The nucleons of the target and projectile interact via two and three-body Skyrme forces 
and Yukawa potential. The isospin degrees of freedom is treated explicitly by employing 
a symmetry potential and explicit Coulomb forces between protons of the colliding target 
and projectile. This helps in achieving the correct distribution of protons and neutrons 
within the nucleus. 

The hadrons propagate using Hamilton equations of motion: 

dfi _ d < H > dfi _ d < H > 

dt dpi dt dvi 

with 

<#>=<T> + <V r >isthe Hamiltonian. 

i % i j>i 

— * — * — * — * 

x fj (r f , p', t)dfdr'dpdp' . 
The baryon-baryon potential V tj , in the above relation, reads as 



(2) 



(3) 



yii(r' — r\ - V ij +V ij +V ij +V ij 

V V ') — v Skyrme^ V Yukawa ^ V Coul^ V Sym 

+ exp(\ r' - f |//i) + ZjZje 2 
(| r' — f \/fj) | r' — f | 

+ tJ-T^.S^-r,). (4) 

Po 

Where \i = OAfm, t 3 = —6.66MeV and t± = lOOMeV. Here Zj and Zj denote the charges 
of the % th and j th baryon, and T 3 l , T| are their respective T 3 components (i.e. 1/2 for 
protons and -1/2 for neutrons). The Meson potential consists of Coulomb interaction 

only. The parameters \i and ti, , t± are adjusted to the real part of the nucleonic 

optical potential. For the density dependence of the nucleon optical potential, standard 

4 



Skyrme-type parameterizations is employed. The Yukawa term is quite similar to the 



surface energy coefficient used in the calculations of nuclear potential for fusion 



13. The 



binary nucleon-nucleon collisions are included by employing collision term of well known 



15 



181 ] . The binary collisions are 



VUU- Boltzmann-Uehling-Uhlenbeck (BUU) equation 
allowed stochastically, in a similar way as are done in all transport models. During the 
propagation, two nucleons are supposed to suffer a binary collision if the distance between 
their centroids 



r, ; - r,- 



<J—, a tot = v(Vs,type) } (5) 



"type" denotes the ingoing collision partners (N-N, N-5, N-7r...). In addition, Pauli blocking 
(of the final state) of baryons is taken into account by checking the phase space densities 
in the final states. The final phase space fractions Pi and P2 which are already occupied 
by other nucleons, are determined for each of the scattering baryons. The collision is then 
blocked with probability 



block 



l_ ( l_ Pl)( l_P 2 ). (6) 



II. RESULTS AND DISCUSSIONS 



In the present calculations, a simple spatial clusterization algorithm dubbed as the 



minimum spanning tree (MST) method is used to clusterize the phase space 



19| . which is 



generated by IQMD Model. We however, also acknowledge that more microscopic algorithm 



routines are also available in the literature 



131 ] . By using the asymmetric (colliding) nuclei, 



the effect of mass asymmetry can be analyzed without varying the total mass of the system. 
We have fixed {Atot = At + Ap = 152) and varied the asymmetry of the reaction just 
like this: 26 Fe 56 + 44 Ru 96 {rj = 0.2), 24 CV 50 + 44 Ru 102 (77 = 0.3), 20 Ca 40 + 50 Sn 112 (77 = 0.4), 
16 S 32 + 50 Sn 120 (r, = 0.5), 14 ^ 28 + 54 Xe 124 (r, = 0.6), 8 16 + 54 Xe 136 (r, = 0.7). 
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FIG. 1: Mean multiplicity of intermediate mass fragments as a function of Zb ounc i at E c . m . = 250 
MeV/nucleon for soft equation of state. Different lines represent the different asymmetries varying 
from 0.2 to 0.7. Here, the values of ^6 oun <i are recorded at different time steps. 

Due to the repulsive nature of Coulomb interactions, one is not able to know the exact 
nature of asymmetry in the reaction dynamics. To understand the role of asymmetry 
beyond the Coulomb effects, we switch off the Coulomb force in our analysis. Additionaly, 
we keep the center-of-mass energy fixed throughout the analysis. 

In order to study the correlation between the < N IMF > and Z bound , it is necessary to 
understand the time evolution of intermediate mass fragments as well as Zb oun d, which is 
shown in Fig. 1. One learns from this figure that the mean multiplicity of IMF increases first 
with the increase of time and then attains equilibrium at later times. The system having 
least asymmetry gives rise to more IMF's as compared to system having large asymmetry. 
This might be due to the reason that as one move towards the large asymmetries, then size 
of the fragments becomes larger than the size of IMFs and hence decrease in multiplicity of 
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FIG. 2: Mean multiplicity of intermediate mass fragments as a function of asymmetry parameter 
rj at different Zbound values at E c , m , = 250 MeV/nucleon for soft equation of state. The lines are 
fitted with an equation y = mx + c, where, m represents the slope of line. 

IMFs is observed. On the other hand, one can see that there is a continious decrease in 
the value of Z bound with time This is due to the decay of compound nucleus into lighter 
particles (i.e free nucleons, LCPs etc.). It means the system is still in non-equilibrium 
state. Moreover, the highly asymmetric system produces largest Zbound because in such a 
case most of the part goes uninteracted. 

In Fig. 2, we show the variation of mean multiplicity of intermediate mass fragments with 
the asymmetry of the system at different values of Zbound- The < Nimf > decreases with 
the increase in asymmetry of the system. This is true for lighter as well as heavier values 
of Zbound- The lines are fitted with equation y = mx + c where, y = N IMF , x = rj, and m 
is slope of the line. The slope values are -1.7, -2.7. -2.6 corresponding to Z boU nd = 16, 20, 
27, respectively. This indicates that for heavier Zb oun d, production of < Nimf > is more 
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FIG. 3: Variation of Njmf and Z bound with impact parameter at E CmTn . = 250 MeV/nucleon. Last 
panel shows the variation of Njmf with Zbound for different asymmetries varying from 0.2 to 0.7. 
Here, the values are recorded at different impact parameters varying from central to peripheral 



one. 



sensitive with asymmetry of the system as compared to the lighter Z bound values. Moreover, 
maximum IMF's are produced at Z bound = 20, indicating the limit of IMFs for a system 
having A = 152. 

Experimentalists studied many times the Z bound dependence of Njmf for symmetric 
5| as well as asymmetric systems. Following these attempts, the detailed analysis with 
asymmetry of the reaction is performed in Fig. 3, where we have plotted the impact 
parameter dependence of Njmf (top panel), Z bound (medium panel), and finally N IMF 
versus Z hound (bottom panel). Due to the low excitation energy E = 250 MeV/nucleon, 
central collisions generate repulsion in a manner so that the colliding nuclei breakup into 



S 



IMFs, whereas for the peripheral collisions, the size of the fragment is close to the size of the 
reacting nuclei, and therefore, one sees a very few IMFs. Interestingly, a rise and fall can be 
seen for nearly symmetric systems, which indicates the possible existence of various decay 
modes from the evaporation (fission) mode to the multifragmentation mode and then to the 



vaporization mode 20j. Moreover, this behavior is observed to be disappear with increase in 
asymmetry of the reaction. On the other hand, in Fig. 3(b) the Zbound is found to increase 
with impact parameter of the reaction. As impact parameter increases, participant zone 
decreases and spectator zone increases, which will lead to the increase in the production of 
heavier fragments and hence increase in the value of Zbound with impact parameter. The 
symmetric systems are more sensitive with the impact parameter dependence of Zbound as 
compared to the asymmetric systems. The maximum asymmetry means that from projectile 
or target, one is heaviest one and other is lightest one. This leads to the possibility of IMF's 
production even at central collisions. The change in geometry does not alter too much the 
production of IMF's in asymmetric systems as compared to symmetric systems. This is 
further elaborated in Fig. 3(c). which shows correlation between the N IMF and Z bound by 
taking into account asymmetry of the reaction (rj = 0.2 to 0.7). Here, the values of Zbound 
are recorded at different impact parameters varying from central to peripheral one. It was 
shown 2j that Zbound allows a very good determination of the impact parameters and hence 
different reaction geometries. The smaller values of Zbound correspond to more central 
collisions. Since the IMF's are produced due to target breakup into pieces, therefore, as 
we vary the asymmetry parameter from r\ = 0.2 to 0.7, the target fragmentation increases. 
The maximum number of IMF's are observed in the Zbound range from 15 to 20. This is in 
agreement with the findings shown in Fig. 2. At lowest asymmetry, we get a rise and fall 
in the production of IMF's with Zbound- But as the asymmetry increases, the curve shows a 
steep variation. These findings are supported by the findings of Fig. 3(a). From this, one 
can see that at highest asymmetry, the size of bounded fragment becomes largest. Therefore 
reaction dynamics changes drastically as one moves from low asymmetry to high asymmetry. 



III. CONCLUSION 



We present our recent results on the fragmentation by varying the asymmetry of the 
reaction between 0.2 and 0.7 at an incident energy of 250 MeV/nucleon. For the present 
study, the total mass of the system is kept constant (A TO t = 152) and asymmetry of the 
reaction is defined by (77 = | (At — Ap)/(Ar + Ap) |). The measured distributions are given 
as a function of the total charge of all projectile fragments, Z bound . We see an interesting 
outcome for large asymmetric colliding niclei. Although nearly symmetric nuclei depict a 
well known trend of rising and falling, this trend, however, is completely missing for large 
asymmetric nuclei. In conclusion, experiments are needed to verify this prediction. 
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